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Abstract 

By proposing an Ansatz for the pressure (measured in terms of the bag con- 
stant) of the hadronic gas in equihbrium, we have formulated a rather simple 
phenomenological model (the extended bag model) which allows one to ana- 
lytically investigate the bulk thermodynamic propertities in the vicinity of the 
phase transition. This makes it also possible to take into account the nonper- 
turbative vacuum effects from both sides of the equilibrium condition. As an 
example of our approach, we have calculated the crossover (critical or tran- 
sition) temperature Tc and the critical chemical potential /Xc (as functions of 
the bag constant) from the noninteracting quark-gluon plasma state equation. 
Our results for TdNf = 0) = 241.5 MeV and Tc{Nf = 2) = 160.6 MeV are 
in good agreement with recent lattice data. The extensive densities such as 
the entropy density, specific heat, etc are calculated as well. A general scheme 
how to calculate the latent heat, the critical energy density, etc within the 
extended bag model is also described and it is applied to the two models of 
the hadronic gas phase. 
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I. The quark-gluon plasma (QGP) phase is a necessary step in the evolution of the 
Exited Matter from Big Bang to the present days. Apparently, the only way to study this 
phase of the expansion of the Universe is nuclear (heavy ion) collisions at high energies 
which makes it possible "to recreate conditions akin to the first moments of the Early 
Universe, the Big Bang, in the laboratory" [1]. Because of the confinement phenomenon, 
the nonperturbative vacuum structure must play a very important role in the transition 
from QGP to the formation of the hadronic particles (i. e., hadronization) and vice versa. 
As it was underlined in our papers [2] , any correct model of the nonperturbative effects such 
as quark confinement or dynamical chiral symmetry breaking (DCSB) becomes a model 
of the true QCD ground state (i. e., the nonperturbative vacuum) and the other way 
around. Thus the difference between the perturbative (always normalizable to zero) and the 
nonperturbative vacua appears to be necessarily nonzero and finite to describe the above 
mentioned nonperturbative phenomena at zero temperature. The existence of the finite 
vacuum energy per unit volume - the bag constant [3, 4] - becomes important for a realistic 
calculation of the transition data from the hadronic gas (HG) to QGP phases at nonzero 
temperature as well. 

There are two main approaches to investigate QGP, namely the resummed finite- 
temperature perturbation theory (effective field theory method) [5, 6] and the lattice one 
[7]. The former breaks down after the fifth order in the QCD coupling constant g because 
of the severe infrared divergences in the Braaten-Pisarski-Kapusta (BPK) series in powers 
of (yf™"*"^" In" (yf [5, 6, 8], but it smoothly incorporates the case of nonzero chemical poten- 
tials). The latter is a powerful nonperturbative tool to calculate equations of state for both 
phases. However, up to now, there are no realistic lattice data available for nonzero chemical 
potential(s) (for problems to introduce it on the lattice see, for example, recent paper [9]). 

As was emphasized in Ref. [10], at present the phase transition at finite baryon chemical 
potential(s) can be only studied within the phenomenological models. The most popular 
among them are, of course, the bag- type models [11], which differ from each other by mod- 
elling the equation of state of the hadronic phase [12, 13]. Within the framework of the 
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bag-type models, a phase transition between the HG and QGP phases is constracted via 
the Gibbs criteria for a phase equihbrium. So the phase transition is necessary of the first 
order. This means that the thermodynamic quantities of interest are discontinuous across 
the critical curve. Our main goal here is to propose a new bag-type model, the so-called 
extended bag model. It complements the standard bag model by an Ansatz for the pressure 
(measured in terms of the bag constant) which allows one to correctly take into account the 
nonperturbative vacuum effects from both sides of the equilibrium condition. We have cal- 
culated the critical chemical potential fic and, especially, the crossover (critical or transition) 
temperature Tc, at which deconfinement phase transition can occur and chiral symmetry is 
restored, in terms of the bag constant. Our numerical results for T(. are in good agreement 
with recent lattice data (see below). The extended bag model makes it also possible to 
analytically investigate the bulk thermodynamic propertities in the vicinity of the phase 
transition. 

II. In the bag-type models the equation of state of the QGP phase is usually approximated 
by the thermal perturbation theory as the ideal (noninteracting) gas consisting of gluons 
and massless quarks. It determines the dependence of the QGP thermodynamical quantities 
such as energy density e and pressure P on the thermodynamical variables, temperature T 
and quark chemical potentials Hf. There exist excellent reviews on the physics of the QGP 
(see, for example, Refs. [14, 15]), as well as on the phase transitions in it [16]. The pressure 
(i. e., the thermodynamic potential Jl, apart from the sign) for the noninteracting QGP is 
given as follows [15] 

P=lhBT' + !^,}T^ + ^y,-B, (1) 

where B is the bag constant (see below), while Nf is the number of different quark flavours. 
In what follows we will consider the values Nf — 0, 1, 2 since the inclusion of the strange 
(s) quark requires a special treatment [17]. The value Nf — describes the case of the pure 
gluon plasma. Note also that the state equation (1) is derived by neglecting quark current 
masses. The constant fsB, entering the equation of state (1), 
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fsB^fsB{Nf) = ^{- + -Nf) (2) 

is the Stefan-Boltzmann (SB) constant, which determines the ideal (noninteracting gluons 
and massless quarks) gas hmit. Obviously, for Nf — it equals to the standard SB constant 
of the ideal gluon gas. It is well known that the energy density e of the noninteracting QGP 
can be obtained from the thermodynamic pothential (1) as follows 

e = 3P + 4S, (3) 

i. e. the bag pressure B determines the deviation from the ideal gas relation. Let us make 
a few remarks in advance. Our calculations of Tc and /ic for the noninteracting QGP are 
not based on the bag model state equation (3). The constraint, determining the phase 
transition, will be obtained with the help of an Ansatz which is beyond the standard bag 
model and it is general (see below, part III). We will use numerical values of the bag constant 
which was obtained from a completely different source, namely it was calculated in the zero 
modes enhancement (ZME) model of the true QCD vacuum on account of the instanton 
contributions as well, (see below, part IV). 

III. The Gibbs conditions for the phase equilibrium between the HG and QGP phases 
a,t T — Tc are formulated as follows [11] 

Ph^Pq^Pc] Th^T^^Tc- 3fXf = iic, (4) 

where subscripts h, q and c refer to HG, QGP phases and transition (critical or crossover) 
region, respectively. At the same time, the difference between Cg — at T = Tc can remain 
finite (nonzero) and determines the latent heat (LH), elh- 

Let us formulate our primary assumption (Ansatz) now. The state equation for the 
hadronic phase (the left hand side of the equilibrium condition (4)) is strongly model depen- 
dent [11-13]. However, in any model the pressure Ph at any values of T and /ih, in particular 
Sit T — Tc and — /J'c , c^-n be measured in terms of the above mentioned bag constant, 
i.e. let us put 
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PH{n,f^c) = -^B, (5) 

O/i + -'V/ 

where the so-called parametric functions bh = bh{Tc, fic) and ah = ah{Tc, fic) describe the 
details of the HG phase at the phase boundary. Evidently, they may only depend on the 
set of independent dimensionless variables (by definition) which characterize the HG phase. 
For example, 

Qh = ah{Tc, He) = Clh{Xc, tc, Vc, Zc), 

bh = bh{Tc, He) = bh{xc, tc, Vc, Zc), (6) 

where 

Xc ^7f7, -Tf^, Vc = — , = /^c-Ro, (7) 

If. If. TTl 

and m denotes the hadron mass while Rq denotes radius of the nucleon, so it allows one to 
take into account finite size effects due to hard core repulsion between nucleons (extended 
volume corrections) [18]. These variables are independent and all other possible dimension- 
less variables are obtained by combination of these, for example, RqTc — Zc/xc, m/Tc — 
Xc/Vc, Hc/B^^^ — Xc/tc,etc. Also the set of independent dimensionless variables at the 
phase boundary (7) may be extended in order to treat the HG phase in more sophisti- 
cated way. However, in any case the parametric functions should be symmetric, i. e. 
ah{Tc:Hc) = ah{-Tc, -^c) and bh{Tc,Hc) = h{-Tc, -He)- 

Prom Eqs. (1) and (4-5) at T = Tc and Hf = /Xc/3, one obtains 

fMN,)Ti^^^T^^^,4=-l^B, (8) 

where we introduced a new "physical" (effective) bag constant as follows 

B = {bh + ah + Nf)B, (9) 

and it linearly depends on Nf as it should be at log-loop level (see below, part IV). In 
connection with our Ansatz (5) a few remarks are in order. Its alternative parametrization 



with respect to Nf, namely Ph{Tc, Hc) — {b'h/a'hNf + 1)B leads to the effective bag constant 
as B = (b'l-^ + a'^,Nf + 1)B. The parameteric functions bh and Uh, as well as b'f^ and a^, as was 
mentioned above, may, in principle, arbitrary depend on dimensionless variables (7). If, for 
example, vanishes a,t — then the linear dependence of B on Nf will be spoiled. In 
other words, the choosen parametrization guarantees the linear dependence of B on Nf and 
the alternative one does not. 

Prom now on Tc and Hc will be calculated in terms of B and not that of old B, i. e. a 
definite numerical value will be assigned to B. So we consider B as the physical bag constant, 
while B as unphysical "bare" one. The bag constant is a universal one and it represents the 
complex nonperturbative structure of the QCD true vacuum. Thus the proposed Ansatz, 
allows one to take into account nonperturbative vacuum effects (parametrized in terms of 
B) from both sides of the equilibrium condition (4). However, it still remains dependent 
on the arbitrary parameter Uh- In order to eliminate this dependence, let us normalize the 
thermodynamic potential at the phase boundary in Eq. (8) to the standard SB constant 
(2). This yields 

fsB{Nf) = {Nf + ah)fsBiNf)^fsB{0) at Nf ^ 0, (10) 

and one immediately arrives at = ah{T(., /ic) — 1- This is our normalization condition 
and it leads to good numerical results for Tc and (see below). So the general constraint 
(8) finally becomes uniquely determined, namely 

f,,T: + ^^T!^^,t-^,S. (11) 

and consequently allows one to investigate the bulk thermodynamic quantites in the vicinity 
of a critical point. Let us emphasize the important observation that the numerical values 
of Tc and /ic calculated from the constraint (11) do not depend on how one approximates 
the equation of state of the hadronic phase. In contrast to the standard bag-type models 
(differed from each other by modelling the hadronic phase) , in the extended bag model their 
values depend only on B which incorporates nonperturbative vacuum effects from both sides 
of the Gibbs equilibrium condition (4) as it has been already emphasized above. 
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The QGP energy density and pressure, however, remain undetermined with our Ansatz 
(5) at this stage. In terms of B they become 



(12) 



and because of Eq. (3), 

^«W) = ^c-(jV, + l)l(iV, + l) + 6/- ('3' 

As mentioned above, the unknown bh reflects the fact that the hadronic phase state equation 
is strongly model dependent. The unknown bh is the price paied to determine the above 
mentioned physical quantities with our Ansatz. Precisely for this reason, the bag model 
state equation (3) plays no role in our numerical investigation of the phase transition with 
the constraint (11) from which Tc, as well as Hc, can be derived. Below (see part VI) a 
general scheme, how to calculate bh within our model, will be developed. Concluding this 
part, let us note that our model cetainly requires the coexistence regime between QGP and 
HG phases since eq(Tc) as given by Eq. (13) explicitly depends on bh which describes details 
of the HG phase at T = Tc. 

IV. Let us discuss the possible values of the bag constant itself now. The bag constant 
is the difference between the energy density of the perturbative and the nonperturbative 
QCD vacua (at zero temperature). The former one can always be normalized to zero, so 
the bag constant is defined as B — —e, where e is the energy density of the nonperturbative 
vacuum. This is always negative, consequently the bag constant is always positive. Let us 
start from the so-called standard value. In the random instanton liquid model (RILM) [19] 
of the QCD vacuum, for a dilute ensemble, one has ej — — (1/4)(11 — |-/V/) x 1.0 fm~^. 
Then the standard value of the bag constant is 

B,t = Bi = -ei = (0.00417 - iV/0.00025) GeV^. (14) 

Thus one can conclude in that the standard value of the bag constant is determined by the 
instanton component of the nonperturbative QCD vacuum only. Note that for Nf — 3 it 



coincides with the estimate of the QCD sum rules approach on account of the phenomenolog- 
ical value of the gluon condensate [20]. From (14) it also follows that the difference between 
values of the standard bag constant for different values of Nf is rather small but not neg- 
ligible. But the main problem with the value of the bag constant, as given by expression 
(14), is, of course, its wrong dependence on Nf. It decreases with increasing Nf that defies 
a physical interpretation of the bag constant as the energy per unit volume. That is why 
the value of the bag constant at the expense of the instanton contributions only is at least 
not complete. 

Instantons are classical solutions to the dynamical equations of motion of nonabelian 
gluon fields which contribute to the nonperturbative vacuum energy density. Therefore 
they are unable to explain confinement phenomenon, which, no doubt, is a quantum non- 
perturbative effect. The QCD vacuum has a much more remarkable (richer) topological 
structure than instantons alone can provide. The dynamical mechanisms of such nonper- 
turbative effects as quark confinement and DCSB are closely related to the complicated 
topological large scale structure of the QCD true vacuum. Assuming that the low-frequency 
modes of the Yang-Mills fields can be enhanced due to the possible nonperturbative IR 
divergences in the true QCD vacuum [21], we have recently proposed the zero modes en- 
hancement (ZME) model of the true QCD vacuum [2]. This is based on the solution to 
the Schwinger-Dyson (SD) equation for the quark propagator in the infrared domain. We 
have shown that this model reveals several desirable and promising features. A single quark 
(heavy or light) is always off mass-shell, i.e. the quark propagator has no poles. It also 
implies DCSB at the fundamental quark level, i.e. a chiral symmetry preserving solution 
is forbidden and a chiral symmetry violating solution is required. We have calculated con- 
tributions to the vacuum energy density at log- loop level [2], coming from the confining 
quarks with dynamically generated masses, €g, and of the nonperturbative gluons, eg, due 
to the enhancement of zero modes. These contributions are almost the same and the sum 
ise = eq + €g = -0.0015 GeV^ - 0.0016 GeV^ = -0.0031 GeV^. Thus the value of the bag 
constant, as given by the ZME model, is 

8 



BzME = -(eg + Nfeq) = (0.0016 + A^/0.0015) GeV"^, (15) 

where we introduced the dependence on Nf since Cq itself gives a single confining quark 
contribution to the vacuum energy density. However, neither contributions (15) nor (14) 
are complete. In the above mentioned papers [2], it was already explained in detail why the 
instanton-type fluctuations are needed for the ZME model. It has been proposed there to 
add Bi, given by Eq. (14), to ZME model value (15) in order to get a more reahstic value 
of the bag constant. Thus one obtains 

B = Bi + BzME = -et = -{ei + €g + Nf€q) = (0.00577 + iV/0.00125) GeV^. (16) 

Numerically our values are as follows 

B{Nf = 0) ~ 0.006 GeV* ~ (278 MeV)* ~ 0.78 GeV/fm\ 
B(Nf = 1) ~ 0.007 GeV^ ~ (290 MeV)"^ ~ 0.91 GeV/fm^, 
B{Nf = 2) ~ 0.008 GeV^ ~ (300 MeV)* ~ 1.04 GeV/fm\ (17) 

We will use these values of the bag constant. Evidently, our value (16) overestimate the 
MIT bag model value of the bag constant [9] at least by one order of magnitude. 

All values of the bag constant below the so-called standard value (14) as well as the 
standard value itself should be ruled out since it does not account for all components of 
the true QCD vacuum, as we pointed out above (see also Refs. [2]). There exist already 
phenomenological estimates [22] as well as lattice calculations [23] preffering a bigger-than- 
standard value of the bag constant. The above described components produce the main 
(leading) contribution to the vacuum energy density and consequently to the bag constant. 
The next-to-leading contributions (as given by the effective potential for composite operators 
at two-loop level [24]) are /i^-order, where h is the Plank constant. Thus they are suppressed 
at least by one order of magnitude in comparison with our values (16). It has been noticed 
in Ref. [1] that nobody knows yet how big the bag constant might be, but generally it is 
thought that it is about 1 GeV/ fw?. The proposed value (17) for the most realistic case of 
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the two thermal quark species {Nf — 2), which is more or less realized in heavy ion collisions 
at high energies, is in agreement with this expectation. 

V. Phase diagram (Tc, /^c) for the physically relevant case Nf = 2, as determined from 
the critical curve (11), is shown in Fig. 1. For end points of this curve (Tc, 0) and (0, //c) 
our data are 

T^{Nf = 0) = 241.5 MeV, 
T^{Nf = 1) = 186.8 MeV, 

%{Nf = 2) = 160.6 MeV, (18) 

and 

li^{Nf = 1) = 1833.8 MeV, 

li^{Nf = 2) = 1441.4 MeV, (19) 

respectively. Thus our value of Tc for the physically relevant case of the two quark species 
A^/ = 2 at /ic = is in fair agreement with its best numerical estimate which comes from the 
observed spectrum of low-pT secondaries in high-energy hadron-hadron collisions, according 
to the analysis of Hagedorn (see Ref. [25] and references therein). 

Let us compare our results for Tc with the finite temperature lattice QCD simulations 
with two light staggered (Kogut-Susskind) quarks represented in Ref. [26]. This can be 
done by means of our data shown in Eq. (18) at /ic = only since, to our best knowledge, 
until now there are no realistic lattice data available for nonzero chemical potentials. The 
lattice result is Tc{Nf = 2) = 155(9) MeV with a systematic uncertanty of about 15%, so 
the agreement with our value (18) is rather good. From our data (18) it also follows that the 
critical temperature Tc for the pure SU{3) gauge theory {Nf — 0) is much higher than for 
the full QCD with two flavors {Nf = 2). This is in agreement with lattice calculations [26, 
27], as well as with arguments of a simple percolation model [27]. A rather small discrepancy 
between our value of the critical temperature for the pure SU (3) gauge theory displayed in 
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Eqs. (18) for Nf ^ and recent lattice calculation Tiatt{Nf = 0) ~ 260 MeV [28] can be 
explained as follows. Prom the constraint (11) at //c = one obtains 



In lattice approach Tc is usually calculated from the string tension in accordance with 
Tj^/a = 0.629(3) [29]. This ratio depends on Nf, i. e. it is smaller for Nf = 2 than 
for Nf — 0, while the string tension itself does not depend on Nf [27]. Using its standard 
value y/a — 420 MeV, one obtains the above mentioned estimate. In a similar way, let us 
evaluate the coefficient in our expression (20) at Nf = 0, while for the bag constant let us 
abstract from its dependence on Nf and use its value for the physically relevant (at nonzero 
temperature) case of the two hght quark species, i. e. B^^'^ — 300 MeV (see Eqs. (17)). 
One immediately obtains T^. — 260.6 MeV in fair agreement with the above mentioned re- 
cent lattice calculations [28, 29]. However, the numerical results of quenched (pure gauge 
theory for staggered quarks) finite-temperature lattice QCD should be (perhaps shghtly) 
reconsidered in the light of the "hard chiral logarithms" problem (see recent review [30] and 
references therein). 

VI. One of the attractive features of our approach is that the thermodynamic quantities 
which are defined as derivatives of the thermodynamic potential Q — —P (extensive den- 
sities, such as the entropy density, specific heat, etc) are uniquely determined at the phase 
boundary. However, let us begin with introducing the metric which shows the existence of 
nontrivial fluctuations in QGP [31, 32]. This is done in two steps. First, the specific heat 
matrix is defined as follows 



The corresponding entropic potential P/T possesses a statistical meaning [31, 32]. The ma- 
trix of its second derivatives Qik — SP'{P/T)/dx'^dx^ with respect to the canonical coordinates 
X* = (1/T, —fj,i/T) defines the avarage scales of the fiuctuations by < gikSx'^dx^ >= 1. There 
is an intimate connection between g^k and the specific heat determinant Ajj. which shows 




(20) 



A 



d^P/dT^ d'^P/dTdiif 

d^P/dTdfif d^P/dfj,} 
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that if the former is positive definite then the latter is also positive definite and the other 
way around [31, 32]. Then the fluctuations remain flnite. For the noninteracting medium 
the specific heat determinant is always positive definite. Using Eqs. (1-2), one obtains 
A{T,Hf) = 4/55^4 + ^N]iJ,j + (6.4 + 1.2Nf)NfT^nj. At the transition phase it finally 



. For 



becomes (on account of Eq. (11)) A{T^, = 4Nf[j^B + (0.177 - 0.133N f)T^ijl 
Nf = 0,1 it is automatically positive definite and, using our numerical results for Nf = 2, 
it is also positive. Thus there are no nontrivial fluctuations in the noninteracting QGP in 
the vicinity of the phase transition as it should be indeed. 

The entropy density is defined as s = , so, on account of Eqs. (1-2), it is, 

s — ^fssT^ + NffjLjT. At the phase boundary it finally becomes 



s{Tc, Pic) = psBT! + ^iATc (21) 

Its behaviour along the critical curve (11) is shown in Figs. 2. 

The specific heat is defined as follows c = T(|^j , and again using Eqs. (1-2), one 
obtains c = ifssT^ + Nf/i'jT. At the phase transition it finally becomes 

c(T„ /.,) = 4fsBT^ + -^filn. (22) 

Its behaviour along the critical curve (11) is shown in Figs. 3. 

The net quark number density is Uf — (^)y- Again using Eqs. (1-2), one obtains 
Uf — Nffj,f{T^ + ^/i/). Since the baryon (B) number of a quark is 1/3, the net baryon 
number density ub = ^Uf at the phase transition finally becomes 

nB{T^,l^c) = ^{T^ + ^yc). (23) 

Its behaviour along the critical curve (11) is shown in Figs. 4. 

One of the important observables measured in heavy ion collisions is the specific entropy 
per baryon, s/ub (see, Refs. [10, 33] and references therein). Its behaviour across and 
along the phase transition may shed light on the strangeness production as a signal of QGP 
formation in heavy ion collisions. At the phase transition it is determined as follows 
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where s{Tc, Hc) and usiTc, fic) are given by Eqs. (21) and (23), respectively. Its numerical 
values can be easily obtained from curves in Figs. 5. In Fig. 1, the path of constant 
s/ub — 50, which is expected from the QGP fireball, is additionally shown. From Fig. 1 it 
also follows that for the QGP fireball temperature T ~ 215 ± 10 MeV the baryochemical 
potential is hb — 340 ± 20 MeV as it should be indeed [33]. 

Let us now compute the numerical values of the entropy density s and specific heat c 
at /Xc = since the dependence on temperature is more important than the dependence on 
chemical potential. From Eq. (21) and constraint equation (11) one obtains 

Numerically this gives 

s^{Nf = 0) = 0.0992 Ge1/^ 
s^(Nf = 1) = 0.0751 GeV\ 

Sc{Nf = 2) = 0.0666 GeV\ (26) 

The numerical values for the specific heat at the phase boundary Cc may be obtained from 
Eq. (26) in accordance with the relation Cc = 3sc, which comes from Eqs. (21) and (22). 

VII. Here let us develop a general method of calculating the parametric function bh{Tc, ji^ 
which remains still unknown in Eqs. (12-13). This makes it also possible to calculate the 
bulk thermodynamic quantities such as the critical energy density, ec, the latent heat, elt-i 
etc, in the vicinity of the transition point T — T^. 

Since the parametric function bh{Tc, f^c) describes details of the HG phase structure, 
in order to determine it, one, evidently, needs to choose the concrete equation of state of 
the HG phase. In this part of our work in what follows, Ph{Tc) will denote the concrete 
equation of state of the HG phase at T = Tc for point-like particles. In order to take into 
account extended volume corrections [18] it should be divided by factor [1 -|- Vonh{T, /i)], 
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where Vq — (47rit!Q/3) with Rq ~ 0.8 /m, is the volume of nucleon, and nh{T,ii) is the 
baryon number density for the point-hke particles, i. e. it should be calculated from P/^ as 
follows 



The more realistic hadron pressure thus becomes 

Ph 



H 



[i + v^on„(r,//)]' 

Solving now Eq. (12) against hh with substitution Ph — > P^, one obtains 



(28) 



l^hM,,.)- + . (29) 

-^^^ B-{Nf + l)PH{%,^ic) ^ ' 

So this relation allows one to calculate hh{Tc, Hc) on account of the chosen equation for the HG 
phase, Ph, by taking the values of Tc and fic from the general constraint (11). Substituting, 
such calculated value of bh{Tc,iJ,c), back to Eqs. (12-13), one obtains the numerical values 
of Pq{Tc) and the critical energy density ec in our model. In order to calculate the latent 
heat, it is necessary first to calculate the energy density of the HG phase as follows 

--^(t).-Kt).-^^- (30) 
Similar to Ph, the more realistic energy density becomes 

^ [l + HMT,/.)]- ^^^^ 
Then the more realistic value of the latent heat is elh — — ^h, where ec is given by Eq. 
(13). This is a general scheme how to calculate clh, Cc, etc within the extended bag model. 
Let us note that our Ansatz (5-7) automatically incorporates extended volume corrections. 

VIII. In this part let us explicitly show how the above formulated general method works 
by approximating the hadronic phase by the ideal (noninter acting) gas of massless mesons. 
The state equation in this case is [11, 12, 34] 

Ph = ^6, = gh^T\ (32) 
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where gh is the number of hadronic degrees of freedom. Evidently, in this case there are no 
extended volume corrections, i. e. Ph — Pr and = en- Prom Eq. (29), on account of Eq. 
(11) at /Xc = and Eq. (32) at T = Tc, one finally obtains 

. (33) 

i. e. in this simple case, the parametric function becomes constant, not depending on T^. 
The critical energy density (13) then becomes 

The meson gas energy density at T = Tc, on account of Eqs. (32) and (11), is 
and the LH becomes 

30/sB - ^hvr^ ~ 

eLH = Cc - = B. 36) 

7.5(A'/ + l)tsB 

Thus these expressions allow one to calculate the bulk thermodynamic quantities in the 
vicinity of the phase transition in terms of the fundamental quantity, the bag constant 
B. The numerical results in the massless pion gas limit are {g^ = 3, Nf = 2): Cc — 
1.358 GeV/fm^, = 0.084 GeV/fm^, elh = 1-274 GeV/fm^, in fair agreement with 
Ref . [34] . This can be explained by the fact that the authors of the above mentioned paper 
a priori use the value of Tc which coincides with the value calculated with our method and 
shown in Eq. (18) for Nf — 2. By construction, the phase transition is of first order, i. e. 
elh is discontinuous (finite). 

Let us now calculate Sh and c^, in this model in order to estimate the jumps in these 
quantities between the HG and QGP phases. The entropy density in the hadronic phase 
is defined as Sh — dPh/dT, where Ph is given by Eq. (32). On account of the constraint 
equation (11) at /Xc = 0, one finally obtains 

(rr\ 9h7r^f 3B ^3/4 

=22:5 ((AT, + 1)/ J ■ (3^) 
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Its numerical value at Qh — ^, Nf — 2 is: Sh — 0.0054 GeV^. The specific heat in the hadronic 
phase is defined as Ch = T(^d'^Ph/dT'^^ . Again on account of the constraint equation (11) at 
lie = 0, one obtains ChiTc) = 3sh{Tc), where Sfe(Tc) is given by Eq. (37). So at gi/j = 3, A^/ = 2 
numerically it is Sh — 0.0162 GeV^. In the HG phase these numbers thus are by one order 
of magnitude less than the corresponding numbers in the QGP phase (see Eqs. (26)). 

It is instructive to compare the numerical value of Tc which follows from the standard 
bag model with that of the extended bag model given in Eqs. (18) for Nf — 2. In the 
standard bag model the pressure, as given by Eq. (1) at /// = 0, should be directly equated 
to Ph, Eq. (32), at T = Tc. This gives the constraint equation as follows 

Tc=(^^ (38) 

Assigning now our value of the bag constant (17) to B, for the massless pion gas limit 
{gh — 3, Nf — 2) one obtains Tc ~ 216 MeV. Of course, this value substantially contradicts 
our value (18) as well as the above presented lattice result (part V). 

IX. Let us consider now a more sophisticated model of the HG phase which consists of 
massless pion gas and nucleons, antinucleons with masses m. The thermodynamic potential 
in this case is [10] 

where K2 is the modified Bessel function of the second kind and the first term describes the 
massless pion gas limit. This is the expression for the point-like particles and Ph on account 
of the extended volume corrections should be obtained from Eq. (28). Using now definition 
(27), one has 

MT,f^) = ^TE^^i^2(/^)sinh(Z^). (40) 
The energy density e^, on account of Eqs. (30) and (39), becomes 

= -T^ + -^T^ E ^-^{m-) - lj;K,{l-)] cosh(Z^), (41) 
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where prime denotes the differentiation with respect to the argumentum. Similarly to the 
previous case, eh is to be obtained from Eq. (31) on account of this expression. These 
relations (39-41) are convinient enough to describe the hot HG and low density matter 
(// — > 0). For cold (T = 0) and baryon dense matter, it becomes a fully degenerate Fermi 
gas at zero temperature. So the finite temperature corrections can be found by using a 
power expansion around T = [10]. 

Expression (29) for hh, in general, becomes 

B[\ + V„n^(T,,ti,}\-(N,+ \)PdT,,^,y ' 
where Ph{Tc,iJi,c) and nh{Tc,ijLc) are given by Eqs. (39) and (40), respectively with the 
substitutions T ^ Tc and n /ic- Thus, in comparison with the previous case (part VII), 
now the parametric function is not simply a constant but it depends crucially on and 
/ic which should be taken from the general constraint (11) in order to numerically calculate 
bh{Tc, He) from this expression. 

Let us calculate, in this model, Cc and e^// at /^c = which allows us to compare numerical 
results with those of the previous part as well as with lattice results (see below) . In this limit 
extended volume corrections disappear and one can use classical statistics for nucleons, i. e. 
only the first term in the series expansion (39) [10]. On account of the constraint equation 
(11) at He — 0, then from Eq. (42), one finally obtains {Nf — 2 everywhere below) 

where rric — m/Tc with m — {mp + raj^^j^ and Tc — 160.6 MeV (see Eqs. (18)). In a 
similar way, from Eq. (41) becomes 

'"^^^^ ^ 7^(l^ ^ ^[^^(^c) - ^X(^c)])^. (44) 

Numerically one get eh — 0.1 GeV/fm^. The critical energy density is to be calculated 
from Eq. (13) on account of the numerical value of bh, obtained from expression (43). This 
is €c — 1.35 GeV/ fm^, so the numerical value of the latent heat is elh — 1-25 GeV/ fm^. 

17 



The value of the bag constant as given in Eqs. (17) for Nf — 2 was usedi n all numerical 
calculations above. These values differ only shghtly from those obtained in the previous 
part. Despite the simple models considered above, these numbers look quite reasonable and 
they may be compared with lattice result especially for Cc, discussed for example, in Ref. 
[35]. The latent heat in this model also remains discontinuous, thereby confirming the first 
order nature of the phase transition. The entropy density Sh as well as specific heat Ch, in 
this model also differ very slightly from those calculated in the previous model. 

An extension of our model to the case of running coupling constant in order to treat 
interacting QGP (which has been already rather tentatively discussed in Ref. [36]) along 
with an extrapolation outside the phase boundary in order to constract equilibrium phase 
transition [10] are subjects of the subsequent papers. 

The authors would like to thank J. Zimanyi, K. Toth, Gy. Kluge, T. Biro, P. Levai, and 
T. Csorgo for useful discussions, remarks and support. 
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FIGURES 

FIG. 1. The phase diagram in the plane (Tcfic) measured in inits of MeV. Here and in all 
diagrams below the curves arc only shown for the physically relevant case of the two light quarks 
Nf = 2 (solid lines). Also shown the path of constant specific entropy per baryon s/ub — 50 
(dashed-dotted line). 

FIG. 2. The entropy density as a function of Tc (upper figure) and fic (lower figure). 

FIG. 3. The specific heat as a function of Tc (upper figure) and /itc (lower figure). 

FIG. 4. The baryon number density as a function of Tc (upper figure) and /Xc (lower figure). 
FIG. 5. The specific entropy per baryon s/ub as a function of Tc (upper figure) and fic (lower 
figure) . 
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